Full alignment of colloidal objects by programmed forcing 
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By analysis and simulation we demonstrate two methods for achieving complete orientational 
alignment of a set of identical, asymmetric colloidal objects dispersed randomly in a fluid. Sedi- 
mentation or electrophoresis in a constant field can lead to partial alignment, in which the objects 
rotate about a common body axis, but the phases of rotation for these objects are random. We 
show that this phase disorder can be removed by two forms of programmed forcing. First, simply 
alternating the forcing between two directions reduces the statistical entropy of the orientation arbi- 
trarily. Second, addition of a small rotating component to the applied field in analogy to magnetic 
resonance can lead to phase locking of the objects' orientation. We identify conditions for alignment 
of a broad class of generic objects and discuss practical limitations. 
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Whenever an object responds to external forcing by 
periodic motion, this response gives potential utility in 
characterizing the object, manipulating it, or using it as 
a probe. A prime example is the response of a nuclear 
spin to a static magnetic field. The response to the field 
enables one to characterize the molecular environment 
of the spin[ ] and to control the quantum state of the 
molecules [2]. Analogs using electrical or mechanical os- 
cillations are well known [3-5]. A further feature of mag- 
netic resonance is that the response of many identical 
spins can be manipulated to be coherent the spin states 
of all the atoms are identical so that they all oscillate 
in concert [6]. Coherence brings further benefits: it en- 
hances the observed signal and it provides further infor- 
mation about the constituent objects [ ]. Moreover, once 
coherence is achieved, the objects can respond as one to 
further external forcing [8]. 

These well-known examples rely on a sharply resonant 
response by the individual objects. Here we examine 
analogous phenomena in the complementary domain of 
purely dissipative responses of identical, asymmetric col- 
loidal objects dispersed in a fluid. In the simple case of 
sedimentation, the object is pulled by gravity acting at 
the center of buoyancy and by hydrodynamic drag forces 
over its surface. For asymmetric objects the drag forces 
generally produce rotation as well as translation. This 
rotational sedimentation effect has aroused recent inter- 
est as a means of organizing colloidal objects [9-1 3]. For 
many objects, the sedimenting force leads to uniform ro- 
tation about a specific direction in the object which aligns 
with the force. Conditions for this "axial alignment" have 
been identified [ ] and related to object shapes [9, 12]. 
Any external field producing motion in the fluid, such as 
electrophoresis, produces analogous rotational effects [14], 
as described below. 

This axial alignment under constant forcing is neces- 
sarily incomplete. Even when all the objects are rotating 
at the same rate about the same body axis, they all have 



arbitrary angular orientations about the aligning axis. 
Here we show how programmed, time-dependent forcing 
can remove this disorder so that all the objects are ro- 
tating together coherently. 

We explain these effects in the simple context of sedi- 
mentation of rigid, asymmetric, noninteracting colloidal 
objects. We first state the equation of motion that gov- 
erns rotation of each object. Then we demonstrate align- 
ment in the simplest case where the force simply alter- 
nates between two different directions. Under mild condi- 
tions this alternating forcing leads to continually improv- 
ing alignment. We quantify this alignment in terms of 
statistical entropy, showing that on average it decreases 
indefinitely. We next demonstrate alignment via a rotat- 
ing transverse force. Finally we discuss generalizations, 
experimental implementations and practical limitations 
of this method. 

For simplicity we consider a set of identical objects sub- 
jected to a sedimentation force F(t) acting at a "forcing 
point" P in the object. We consider the regime of creep- 
ing flow [ ], in which inertial forces are negligible and 
the force transmitted to a moving object by the medium 
is proportional to the object's velocity. A rotating rigid 
body with center of mass velocity v and angular velocity 
Q experiences a proportional hydrodynamic force F and 
torque r. In terms of the hydrodynamic radius R of the 
object, the proportionality can be written in terms of a 
dimensionless block matrix: 
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where n is the viscosity of the fluid. For simplicity below 
we will use units such that 6ttt] and R are unity. By 
choosing the forcing point P as our origin, we remove 
any external torque, so that any rotation arises entirely 
from the 3x3 "twist matrix" T: Co = TF. The change of 
T owing to this uj can be written [9]: 

T=[uj x ,T], (2) 
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Here we use the notation uj x to denote the antisymmetric 
tensor corresponding to the vector 0: [uj x ]ij = —€ijk ook- 
Alternatively, in a rotating frame that is fixed in the ob- 
ject, T becomes constant and F rotates: F = —uo x F. 

In what follows we restrict our attention to a sub- 
class of T's — denoted "axially aligning" — that align in a 
unique direction, independent of initial orientation, un- 
der constant forcing. Such T's have sufficiently large an- 
tisymmetric parts that they have only one real eigenvalue 
A3 and eigenvector #3. Then T aligns with its #3 along 
F and rotates at angular velocity A3F [16]. An example 
of one such object, made by four conjoined spheres, is 
shown in Fig. la [17]. 

We first consider a minimal forcing program: a mere 
switch of the forcing direction by a "rocking angle" 6 
from its initial direction along the z axis. Fig. lb il- 
lustrates the result of this switch. An initially uniform 
distribution of phase angles <p becomes nonuniform; the 
transient relaxation has reduced the disorder. 

More explicitly, we may define a basis of unit vectors 
ei, £2, £3 in the object whose £3 axis is the direction that 
aligns with the force F, initially along the z axis in the 
lab frame (Fig. la). At a given moment this object's 62 
vector makes an azimuthal angle <p with the lab's y axis. 
This (j) increases steadily in time as explained above. The 
0's of the different objects are presumed to be uniformly 
distributed. After F has switched into the x-z plane and 
the objects have re-aligned, their £2 vectors again differ 
only in their azimuthal angle with the y axis, which we 
denote ip. Evidently an object's angle ^ at a time t after 
the switch depends on the angle <p immediately before 
the switch: ip = ipt{<j>)- As the figure illustrates, the ip 
angles are no longer uniformly distributed: some ip's have 
bunched closer together; others have spread apart. For 
convenience we shall choose a time t so that ^(0) = 0. 
Evidently if the rocking angle 8 = 0, the motion is a 
continuation of the uniform rotation without any switch: 
ip{(p) = (p. If increases from zero by a small amount, 
ip{(p) remains close to (p and ip{(p) remains monotonic [19]. 
Further, for general 0, as <p advances by 27r, ip must ad- 
vance by the same net amount. Thus 2tt = § d(p {dip/dcp) 
[19]. Typical ip{4>Ys are shown in Fig. lc-e. 

We may quantify the bunching effect of ip using the 
probability distribution function p{(p) measured after 
the objects have aligned. The net effect on the distri- 
bution can be quantified using the statistical entropy 
H = — J plrip. This H is maximal for uniform proba- 
bilities and is small when the probability is concentrated 
into small regions [20]. After one switching process, each 
angle <p evolves into some <p. The passage from <p to <p 
involves two steps: we first wait a random fraction of the 
rotational period, so that <p undergoes a shift by a ran- 
dom angle a. We then switch the force F and allow the 
objects to realign. The resulting angle <p = ip{(p + a). 
The corresponding probability distribution p{<p) is found 
using p{(p) dcp = p{(p) d(p, so that p{(p) = p{<p) /ip' {<p + a) . 
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FIG. 1. a) The conjoined-sphere object referenced in the 
main text axially aligned to a downward force F, viewed from 
slightly below the horizontal. A lab basis x-y-z and a body 
basis ei, £2, £3 , centered at the forcing point, are shown. The 
azimuthal angle cp is 0. The object's T matrix in ei,e2,e3 
co-ordinates was computed using the HYDROSUB software 
package [ ]. A colored dot marks the £2 axis, b) Motion 
after a switch in force direction. The red line illustrates the 
trajectory of the colored dot after the forcing direction has 
switched to the x direction. Eleven other trajectories for ob- 
jects initially rotated by different 0's are also shown. After re- 
alignment the points have become non-uniformly distributed 
on the new circle, c), d) e): sample results of repeated rock- 
ing by 90 degrees as described in the text. Left-hand graphs 
show the mapping function ip{<j>). Center graphs show the 
angles <p for 100 randomly oriented copies of the object after 
100 switches. Right-hand graphs show the progression of the 
entropy H with repeated rocking. Final low values are consis- 
tent with the noise floor of the numerics . c) shows the results 
of a monotonic ip function. Line indicates the predicted slope 
from (4). d) shows the results for the object shown in a, e) 
shows the results for a strongly non-monotonic ip{(j>) whose 
entropy does not decrease. 



The change of H resulting from this process can readily 
be computed provided ip' {(p) is positive for all <p[ ]. The 
change AH depends on the shift angle a: 



p(0)ln(^(0 + a)) #. 



(3) 



We now consider the net change of H after a sequence 
of many switches of F. Then the a-averaged change of 
H, (AH) can be written 



(AH)= J>ip(4> + a)) In 



(4) 



The (p{<p + a)) is an unimportant positive constant. 
The remaining integral is a constant that is necessarily 
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negative [20], owing to the convexity of the logarithm. 
Thus after many switches, the entropy must decrease in- 
definitely on average and the probability p becomes con- 
fined to arbitrarily small regions of <j). As Fig. lc and d 
illustrate, the initial set of 0's often evolve into a single 
final (f) [21]. Further, this reduction of H often occurs 
even when ip' is not positive definite (Fig. Id). Given 
a mixture of two or more alignable species, this method 
aligns each species. However, one readily finds examples 
"0(0) where the entropy does not decrease (Fig. le). In 
any case the direction of alignment is not controlled in 
this method. 

A second method of forcing can achieve alignment in a 
controlled direction. We add a rotating transverse force 
to the original static force, so that the force vector is 
tilted at an angle from the z axis. Thus F rotates at a 
constant angular velocity SI = \Sl\ z. 

F(t) = \F\ ( z cos 6 + [x cos Sit + y sin Sit] sin ) . (5) 

In a co-rotating frame rotating at angular velocity SI, 
the force F becomes constant. With a proper choice of 
SI and the objects too may evolve into a state of co- 
rotation with the force, with a common orientation. Co- 
rotation requires that T remain in an orientation such 
that SI = TF. As noted above, a fixed force with 6 = 
leads to co-rotation with SI = XsF. In the ei, £2, £3 body 
reference frame, denoted with subscript 6, co-rotation 
means Sl b = T b F b for some constant F b and correspond- 
ing Sl b . The forcing parameters required are evidently a) 
\F\ 2 = |F 6 | 2 , b) \Sl\ 2 = \Sl b \ 2 and c) cos6> = F b • Sl b . 

With this same choice of parameters other co-rotating 
F's are also possible. We denote these as F^. Because 
of condition a) such F^s have the same magnitude as F, 
so that the only difference in the F b s is in their direc- 
tion Ffr. In terms of the matrix T5, condition b) can be 
written F 6 ' T 6 T T 6 F 6 ' = F b T b T T b F b ^(= |^| 2 /|F| 2 ). (If 
Fl satisfies this condition, so does —F^.) This condition 
restricts F 6 ' to two (closed) curves on the unit sphere like 
the dashed line in Fig. 2a-c. Likewise condition c) reads 
H T 6 H /\T b Fl\ = F b T b F b /\T b F b \ (= costf). This 
condition restricts F^ to a second pair of curves on the 
unit sphere. Any intersection of these curves represents 
an Fl that co-rotates with the given forcing. 

These compatible F b s can readily be found when the 
tilt angle 6 is small. Figs. 2a-c show the behavior of the 
co-rotating F^ s as one increases 6 from with \Sl\ = X3F. 
The condition-b) curves, enforcing the magnitude of SI, 
are thus independent of 6 and do not change. These 
curves pass through £3 and —£3, i.e., the F b for 6 = 
0. Condition c), enforcing cos# requires F b || £3 when 
6 = 0. As 6 increases, the condition-c) curves expand to 
small rings encircling ±£3. Each ring must intersect its 
condition-b) curve twice. The four co-rotating F^ s are 
then two adjacent pairs near ±£3. In what follows we 



consider only the two intersections adjacent to the stable 
+£3 direction 

To achieve full alignment, an arbitrary initial state 
must evolve into one of these two co-rotating states. We 
argue that this alignment occurs generally for sufficiently 
small 6. We consider the motion of T in the co-rotating 
frame, in which F and SI are fixed. We first align T axi- 
ally using a constant F. We can express the orientation 
of T using a rotation vector 77. We distinguish angular 
rotations 773 along the SI axis from rotations n± perpen- 
dicular to it. In this initial state, rj± = 0, while the 
orientation 773 about the aligning axis is arbitrary: axial 
alignment has restricted ff to a closed one-dimensional 
curve of possible values. A small rotation 773 has no fur- 
ther effect on the motion, but a small rotation n± leads 
to a stable return to n± = 0. 

Now we increase to a small nonzero value. Then ff is 
no longer constant. The time derivative ffo(ff) differs from 
770 (ff) by a small, smooth perturbation. We now suppose 
that n± converges to some 773-dependent value near 0. 
The remaining motion, is along a one-dimensional closed 
curve, slightly distorted from the neutral curve found at 
6 = 0. We may express this residual one-dimensional 
motion using the co-ordinate 773, whose time derivative 
is some function 7730(773). The two fixed points identified 
above are necessarily fixed points 77I, with 7)3(77!) = 0. 
Since fjse is a smooth function on the curve, its deriva- 
tives at these two fixed points are necessarily opposite 
in sign; hence opposite in stability. Thus in this picture 
all 773 must must converge to the stable 773. Numerical 
studies like those of Fig. 1 confirmed this finding for 
numerous asymmetric T's . 

A local stability analysis confirms the opposite stabil- 
ities of the two fixed points F b near e%. Starting from a 
given forcing with a corresponding F b , we rotate T by a 
slight angular displacement ff from its fixed-point state 
T*, in the co-rotating frame so that T = T* + [77 x , T*]. 
The undisplaced state ff = is co-rotating, so that ff =0. 
Near ff = 0, ff varies linearly with 77: 77 = K77 for some 
"stability matrix" K = [T*F X - (T*F) X ] [ ]. From this 
K one may determine by standard methods [22] whether 
the motion returns stably to the co-rotating state at T* 
after the small displacement 77; and, indeed, one finds 
that the two fixed points of the previous paragraph have 
opposite stabilities. Figs. 2a-c use this K to determine 
which F b represent stable co-rotating states. The figure 
shows large stable and unstable regions of F b . We note 
that the aligning axis £3, being neutrally stable, lies at 
the boundary between the stable and unstable regions. 

These alignment methods are applicable for any situa- 
tion where an object responds to a vector field by rotat- 
ing. For example, any asymmetric object with an elec- 
trophoretic mobility has such a rotational response[14]. 
Any proportional response is necessarily governed by a T 
matrix like that above. This is true for deformable ob- 
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FIG. 2. a) The sphere of possible directions of the body frame 
force for the object of Fig. la. The solid curve gives the 
Fl s satisfying constraint c in the text, enforcing an angle 
— 0.3. The dashed curve gives the F's satisfying constraint 
b in the text, enforcing \Q\ = A3 \F\. The intersections of the 
two curves mark co-rotating states. The light shading marks 
regions of stable fixed points Fb for this object, b), Same, for 
the T used for Fig. 1 e with = 0.2. c) Same as b, with Fb 
chosen to lie far from the aligning direction £3 to illustrate 
alignment with large tilt angle. Here \Q / A3 |i^| ; instead, 
fib = ^bFb. d) Time sequence of orientations using the object 
and alignment protocol described in a. Three arbitrarily ori- 
ented copies of the object, represented by their ei, eb, £3 bases 
are shown. Colored arrows indicate the direction of forcing. 
The three objects evolved to a common final orientation, e) 
Similar time sequence for the T matrix used for Fig. 1 e. f) 
Similar time sequence for the T matrix and forcing condition 
shown in c. 



jects whenever the driving force is weak enough to avoid 
deformation. It is true for fluctuating shapes whenever 
the driving is weak enough to average over the fluctua- 
tions of T. 

Complete alignment brings benefits that are not pos- 
sible in axial alignment. Once a set of objects have been 
completely aligned, they must respond identically to sub- 
sequent forcing. Then eg. by gradually increasing the Q 
of Eq. 5 one may achieve phase locking without knowing 
about the T of the objects in advance [ ]. Further exten- 
sions of programmed phoresis hold promise for orienting 
non-axially-aligning objects not considered here. 

In practice, alignment is degraded by rotational diffu- 
sion, characterized by a relaxation time td, which scales 



as the cube of the object's size R [ 5]. The driving forces 
must be sufficiently strong to produce rotational speeds 
uj much larger than 1/td- Thus alignment is strongly 
degraded as R is decreased. The object of Fig. la has 
Qtd ~ 260. Under electrophoresis an object with a nom- 
inal mobility of 10 -8 meters/ (volt sec) and field of 10 4 
volts/m, would need to be over 10 microns in size to give 
comparable Qtd- Thus these methods require objects of 
near-micron scale or larger. 

Unless the objects are very dilute, their mutual hy- 
drodynamic interactions would be significant: an object 
B is advected by the perturbed Oseen flow around a 
nearby object A [ ]. Only the gradient of this flow ro- 
tates object B. The rotational perturbations on an ob- 
ject are thus shorter range and hence weaker than are 
translational hydrodynamic interactions. Ordinary elec- 
trophoresis produces no Oseen flow in the host fluid [14]; 
thus electrophoresis is less subject to interaction effects. 
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In a 10 4 g centrifuge this object with spheres of 1 micron 



